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I. Introduction 

Users in a multihop wireless network convey information to each other with the help of in- 
termediary routing nodes. Generally, packets are forwarded by the relays towards their respective 
destinations in a decode-and-forward fashion. After the seminal paper by Ahlswede et al.[l] on 
network coding for wireline networks, it is known that better performance is possible if intermediate 
nodes are allowed to change the content of their packets [2], [3]. Wireless network coding can 
improve throughput and reliability due to the broadcast nature of the wireless medium, and the 
resulting opportunities to gather information from all audible transmissions [4], [5], [6], [7], [8], 
[9]. Since two-way traffic is inherent to peer-to-peer communication, the two-way relaying channel is 
a key building block for information exchange over multiple hops. This paper provides fundamental 
characterization and limits of the two-way relaying channel with (and without) multihop wireless 
network coding. 

Specifically, we consider two multihop network coding protocols as illustrated in Fig. [TJa) in this 
work, i.e., MAC -Layer network coding (MLNC) and PHY- Layer network coding (PLNC). MLNC 
is a coding operation that happens at or above the media access layer [10], [11], [12], [4], [8], [7]. 
PLNC is based on the coding protocol proposed in [13], [14], [15]. In order to perform MLNC and 
PLNC, the two source nodes respectively transmit their packets to the relay node in the first two 
time slots. Then the relay node constructs a new network-coded packet from the received packets 
and broadcasts it to the two source nodes in the third time slot. If MLNC is adopted - whereby the 
packets are manipulated before channel coding - the transmission rate of the relay node is limited 
by the smaller channel capacity [13]. Since the network coding procedure of PLNC is performed 
instead on channel coded data, PLNC can individually achieve the broadcast channel capacities from 
the relay node to the source/destination nodes [16], [14], [15], [13], [8]. Therefore, PLNC always 
achieves better throughput than MLNC, especially when the two channels are highly asymmetric. 
PLNC's performance gains come at the cost of transceiver complexity. 

A. Motivation and Related Work 

Traditionally (pre-network coding), information exchange between two users via a relay has been 
accomplished by a time-division multihop (TDMH) protocol in four time slots[jl as shown in Fig. 

'Frequency-division could be used as well. In this paper the comparisons focus on time-division systems only, and comparisons 
can be applied to frequency-division systems similarly. 
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[Da). Intuitively, the MLNC and PLNC protocols save one time slot. In [6], two-way relaying for 
cellular systems was considered, while [4] and [17] proposed an MLNC algorithm effective for 
wireless mesh networks in heavy traffic. The network coding protocol in Fig. djb) can be further 
reduced to two slots if advanced joint coding/decoding - i.e. analog network coding (ANC) - is 
allowed. In this case, both source nodes send their packets to the relay node simultaneously during 
the first slot; then the relay node either amplifies and broadcasts the signals, or broadcasts the XOR- 
ed packets after decoding by successive interference cancellation [18], [19], [9], [20]. The achievable 
rates for analog network coding were studied in [21], [9]. 

Although MLNC and analog network coding have been shown to achieve throughput gains in 
certain environments [6], [19], [7], it is unclear how channel realizations influence the achievable 
rates and whether MLNC and/or PLNC are always better than TDMH in terms of end-to-end 
throughput. In terms of implementation complexity, the four techniques can be ranked from low 
to high as TDMH, MLNC, PLNC ANC. Analog network coding requires stringent synchronization 
and joint decoding and could suffer from decoding error propagation due to channel estimation and 
quantization errors. Characterizing the exact achievable rate regions of TDMH, MLNC and PLNC 
will provide guidance on the tradeoffs between them. 

A frequently neglected but very important consideration for network coding is that the traffic 
patterns can vary significantly, and have a significant effect on the ability to achieve the gains 
promised by network coding. For example, data downloads are essentially one-way traffic, while 
peer-to-peer conversations are fairly symmetric. Analyzing the gains from network coding in a two- 
way relay network in view of the traffic pattern allows insight into how such systems should be 
designed. In this paper we explicitly consider the traffic pattern in our results. 

Finally, network coding can be used to exploit cooperative diversity between source and destination 
nodes[22], [23]. Since network coding is able to provide diversity as well as throughput gain, it 
is of interest to understand the diversity-multiplexing tradeoffs (DMTsj^l of MLNC and PLNC and 
determine if they are better than TDMH's. Since we consider two-way transmission over multiple 
relays, this plurality of relays may cooperate in a number of different ways or not at all, and each 
cooperation scenario leads to a different DMT result for TDMH, MLNC and PLNC. 

2 The diversity-multiplexing tradeoff (DMT) for point-to-point multiple input and multiple output (MIMO) channels was found in 
[24], and has become a popular metric for comparing transmission protocols. 
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B. Contributions 

In this paper, we first characterize the exact achievable rate regions of TDMH, MLNC and PLNC 
and show that MLNC is not always superior to TDMH, while PLNC has the largest rate region 
among the three. An opportunistic scheduling algorithm is then presented to achieve larger rate 
region than those achieved independently by TDMH and MLNC. The stability of this scheduling 
under Poisson arrivals is shown via queuing analysis. Subsequently, the sum rates of the three 
protocols with or without traffic constraint are determined. The optimal sum rate with or without 
a traffic pattern constraint - defined as the ratio between the rates of the two directions - is also 
characterized. 

The DMTs of the three transmission protocols with different relay collaboration scenarios are all 
derived. They are quite different from previous multihop DMT results (e.g. see [25], [26], [27], [28]) 
due to their dependence on the traffic pattern, time allocation to each transmission direction, and 
the number of cooperative relay nodes in the system. We first characterize the DMT of TDMH with 
relay collaboration and show that it coincides to that of TDMH with an optimally selected relay 
node. Then, MLNC and PLNC are shown to achieve the same DMT in the high signal-to-noise 
ratio (SNR) regime, and that better tradeoffs are attained if an optimal relay node is selected to 
broadcast. Considering these results, we conclude that selecting an optimal relay node to receive 
and transmit (or broadcast) is preferable. Finally, we find that MLNC and PLNC can have a worse 
DMT if there is suboptimal time allocation for a certain traffic pattern. Intuitively, if the offered 
traffic load is much higher in the forward direction than the backward direction relative to one of the 
source nodes, then network coding may not be helpful for that source since it presumes a symmetric 
data rate in the high SNR regime. 

II. System Model and Assumptions 

Consider a multihop wireless network in which information exchange by multihop routing can 
be characterized by considering the two-way relaying system illustrated in Fig. [T] where the two 
source nodes A and B would like to exchange their data packets Wa and Wb, denoted as two 
binary sequences. In the sequel, we assume there is no direct channel between the two source 
nodes, otherwise, mutihop is not needed [29], [30]. All nodes in the network are assumed to be 
single-antenna and half-duplex, i.e., no nodes can transmit and receive at the same time. Also, we 
assume that channel gains are constant during the whole transmission, and the channel gain between 
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any two nodes X and Y, denoted by h X y, is reciprocal and modeled as a zero-mean, independent, 
circularly symmetric complex Gaussian random variable with variance 1/gxy- 

The core idea of MLNC is that the relay constructs a new network-coded packet Wd after receiving 
Wa and Wb, and then it broadcasts Wd to both source nodes. Wd is obtained by directly XOR-ing 
Wa and W B bitwise, i.e., W D = W A © Wb- As a source node receives W D , it can decode the new 
packet by XOR-ing Wd with its side information, i.e., Wa or Wb- In contrast to XOR-ing data 
contents in MLNC, PLNC does network coding on channel codes, i.e., after channel encoding to 
individual data contents [13], [14], [15]. Taking binary symmetric channel for example, the relay 
broadcasts Xab ® X B a after receiving Wa and W B , where Xab and X B a stand for the channel 
codes of Wa and W B respectively, as shown in Fig. [2] The channel encoder of Wa (Wb) which 
generates X AB (X B a) is designed according to the channel condition from relay node D to source 
node A (B). So when the XOR-ed channel code is received by a source node, it can subtract Xba 
or Xab before channel decoding, as shown in Fig. [2] As shown in [13], [14], [15], as long as 
there exists a common input distribution at the relay node that achieves channel capacities for both 
directions (e.g., binary symmetric channel and AWGN channel), similar idea works. In this paper, 
we assume the existence of such common input distribution. Due to the different constructions, 
the available broadcast rates for MLNC and PLNC are different. For MLNC, the broadcast rate 
is limited by the smaller broadcast capacity, due to the fact that both ends need decode the same 
(XOR-ed) message. For PLNC, each direction can achieve its individual channel capacity [13]. Let 
Cxy denote the channel capacity from nodes X to Y. The broadcast rates of MLNC and PLNC can 
be concluded from [13] [15] as follows. 

Lemma 1: The achievable broadcast rate of MLNC for the relay node is C m i n for both directions, 
where C m - m = mm-fCTo^, Cdb}- If there exists a common input distribution that maximizes the 
mutual information from relay node D to nodes A and B, then the achievable broadcast rates of 
PLNC for both directions are respectively Cda and Cdb- 

In this paper, we also consider the information exchange between two source nodes can be 
completed with multiple relay nodes, as shown in Fig. [EJb). Dab denotes the set of the relay nodes 
available^ between source nodes A and B. Denote by \T>ab\ the number of relay nodes in "Dab, 
which is usually a random variable for different time slots; however, to facilitate the analysis here 
we assume it remains constant during the period of exchanging packets. We assume all nodes in 



3 where "available" means any relay node in Dab can successfully decode the information from both source nodes. 
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Dab are close and able to collaborate under reasonable communication overhead so that every relay 
node can share its received information with other relays. In this context, Dab virtually becomes 
a big relay node equipped with \Dab\ antennas. The channels from node A to Dab become a 
single-input-multiple-output (SIMO) channel (or a MISO channel from Dab to node A). Therefore, 
receive maximum ratio combining (MRC) and transmit MRC can be accomplished in Dab for 
TDMH assuming joint processing can be carried out. Although receive MRC can be performed in 
the first two transmission stages for MLNC and PLNC, it is hard to achieve bidirectional transmit 
MRC in the broadcast stage. Thus in analyzing MLNC and PLNC with relay collaboration we 
consider two scenarios of broadcasting, i.e., all relay nodes broadcast at the same time and only an 
optimally selected relay node broadcasts. 

III. Achievable Rate Region, Opportunistic Network Coding and Scheduling 

In this section, we are interested in determining the end-to-end rate pair (Rab, Rba) achieved 
by the aforementioned three protocols. For convenience, we call Rab the forward rate, Rba the 
backward rate and fx the traffic pattern parameter which is the ratio Rab I Rba- Here we only 
characterize the end-to-end rate regions achieved by TDMH, MLNC and PLNC for the single relay 
network in Fig. Qla) since they are easily extended to the multiple relay case. Given the achievable 
rate regions, two opportunistic packet scheduling algorithms are respectively proposed for MLNC 
and PLNC, and their stability with random arrivals are characterized as well. 

A. Achievable Rate Regions for Two-Way Transmission Protocols over a Signal Relay 

For the two-way relaying system in Fig.[H we assume that Cda and Cdb are achieved by the same 
input distribution. The achievable rate region is basically constructed by the forward and backward 
rate pairs (Rab, Rba)- First consider TDMH that needs four time slots to exchange packets. Since 
its achievable Shannon rate pairs are constrained by time allocations in the four time slots, its 
achievable rate region is 

^tdmh — I (Rab, Rba) '■ Rab < {^iCad, ^Cdb}, Rba < {^Cbd, ^Cda}, = 1 >, (1) 

^ k=i ' 

where {A^} E [0, 1] are the time-allocation parameters for four transmission time slots. Define 

^ab — (1/Cad + ^-/Cdb)^ 1 and S BA = (1/C B d + V^da) -1 - Therefore, we have the following: 
Theorem 1: "Jl TDMH is the triangle with vertices 0, (Sab, 0) and (0, T, B a), as shown in Fig. \3\ 
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Proof: Vertex (0, T> BA ) corresponds to the case of one-way backward traffic, it is achieved by 
setting Ai = A2 = 0, A3 = and A4 = Similarly, vertex (Sab? 0) corresponds to the case of 
one-way forward traffic, achieved by setting A 3 = A4 = 0, Ai = §^ and A 2 = Since CR TDMH 



is 



described by linear constraints, it is convex and thus achievable. Now we show that ft TDMH is also an 
outer bound for TDMH. Consider the four linear constraints of transmission rates in ©. Dividing 
each of them by their corresponding channel capacity and adding them up, we obtain 

Rab Rab Rba Rba _ Rab Rba < \ ^ > . 
^ad ^db <^bd <^da ^ab ^ba ~~ 

k=l 

This is exactly the region below the line connecting vertices (0, T* BA ) and (Sab, 0). ■ 
For MLNC, by Lemma [TJ its achievable rates are defined in similar fashion as in (OQ) as follows 

^mlnc — \ (Rab, Rba) '■ Rab < {AiC^d, A 3 C min }, -R^a < {^Cbd, ^Cmm}, / ] — 1 r • (2) 
^ k=i ' 

Define H AA ± (l/C AD +l/C DA )-\ Y, BB 4 {1/C BD +1/C DB )- 1 and Y, ABB 4 (l/E AB +l/C BD )-\ 

Then we have the following theorem characterizing the achievable rate region of MLNC. 

Theorem 2: If Cdb < Cda, then %. M lnc is the quadrilateral with vertices 0, (J2 AB ,0), (0, Ebb) 

and (Eabb, Eabb), as shown in Fig. \3\a). If Cdb > Cda, then "31 M lnc is the quadrilateral with 

vertices 0, (E AA , 0), (T, ABB ,T, ABB ) and (0, T, BA ), as shown in Fig. \3\b). 

Proof: First consider the case C D b < C D a- The achievable rate region CR M lnc m © becomes 

3^mlnc — \ (Rab, Rba) '■ Rab < {^iC A d, ^Cdb}, Rba < {^Cpp, A 3 (7db}, A& = 1 >. (3) 
^ fc=i J 

Vertex (T, AB ,Q) corresponds to the case of one-way forward traffic, and is achieved by setting 
Ai = A 2 = and A 3 = Vertex (0, Yj BB ) corresponds to the case of one-way backward 
traffic, and is achieved by setting Ai = 0, A 2 = ^P^- and A 3 = Finally vertex (Eabb, Sabb) is 
achieved by setting Ai = , A 2 = ^Cbd an( * ^ 3 = ^Cdb ' ^ ecause @ is a convex quadrilateral 
defined by the three vertices and (0, 0), the region © is achievable by time-sharing among the four 
vertices. Next, by time-sharing we show the quadrilateral CR MLNC in © is also an outer bound for 
MLNC. Consider the four linear constraints in §3$ and divide each of them by their corresponding 
channel capacity. We have 
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Rab , Rab . Rba Rab . Rba \ > , , , , 

77— + TT- + TT- = y + - l^ X k = !> ( 4 «) 

^AD ^DB ^BD ^AB ^BD T~T 

Rba , Rba , Rab Rba . Rab . , 1 , A1 , 

-7T- + TT- + TT~ = v + TT~ < }^ Xk = L {Ab) 

<~>BD ^DB ^AD ^BR ^ AT) 



k=l 



s 



Equation (l4ol) corresponds to the region below the line connecting (E^ B , 0) and (T: abb, ^ abb), 
and (|4FI) corresponds to the region below the line connecting (T^abb, ^abb) and (0, Ebb). This 
completes the proof for the case Cdb < Cda- The proof for the case Cdb > Cda is similar. ■ 
Remark 1: According to Theorems [TJ [2] and Fig.[3ta)(b), MLNC is not always better than TDMH. 
For example in Fig. Oa), MLNC is worse than TDMH when Cdb < Cda and /i is greater than 
certain value. This is because the broadcast rate of MLNC is limited by the worse broadcast channel, 
as stated in Lemma [TJ A hybrid protocol of time sharing between MLNC and TDMH, as indicated 
in Fig. [3ta)(b) by a dashed line, can achieve a larger rate region CovH(3? TDMH , 3^mlnc)> the convex 
hull of Ktdmh and 3? M lnc- Thus, in practice the relay can decide to use TDMH or MLNC according 
to the results in Fig. [3] 

For PLNC, its achievable rate region is constructed as follows according to Lemma [TJ 

^plnc — S (Rab, Rba) '■ Rab < {^iCad, ^Cdb}, Rba < {A 2 Cbd, X^Cda], X k = 1 >. (5) 
^ fc=i J 

Define Y^aba = (1/Eab + Cda/Cbd/CdbY 1 and E BAB 4 (1/E BA + Cdb/Cad/CdaY 1 . Then 
we have the following theorem. 

Theorem 3: l Ji PLNC is the quadrilateral with vertices 0, (£ab,0), {J^aba^bab) and (0, Hba)- 
Proof: The proof is omitted here since it is similar to the proof of Theorem [2l ■ 

Remark 2: Given the above results in Theorems UJ3] and illustrated in Fig. [3l we know the 
achievable rate region of PLNC is larger than CovH(ft TDMH , CR MLNC ) as long as Cda 7^ Cdb- Also, 
it should be noticed that time sharing between TDMH and PLNC does not help achieve larger rate 
region as in the case of MLNC and TDMH. 

B. Opportunistic Network Coding and Scheduling 

The achievable rate regions for the three transmission protocols have been characterized in Section 
IIII- A I In this subsection we investigate the following question: If packets arrive at the two source 
nodes according to a random process, how should the packets be scheduled for transmission to 
maximize the rate region in which the queues are stable? We assume that the buffer size for queuing 
packets at the two source nodes is infinite. Packets of length i (bits) arrive at the queue of source 
node A according to a Poisson process with rate Rab, while packets of length £ (bits) arrive at the 
queue of source node B according to another independent Poisson process with rate Rba- We also 
assume that i is large and channel coding is perfect. 
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We first consider the case where the relay node uses MLNC to route the packets. Consider the 

time right after the previous transmission is complete. Let qA and qs denote the number of packets in 

the queues at node A and B, respectively. We propose the following opportunistic packet scheduling 

algorithm for MLNC. 

Algorithm 1 : Opportunistic MLNC and Scheduling 

Step 1 (qA ■ 1b 7^ 0): Two packets (one from each queue) are sent over Relay D by MLNC. 

Step 2 (qA 0, qs — 0): One packet from node A is sent over Relay D by TDMH. 
Step 3 (qA = 0, qs 7^ 0): One packet from node B is sent over Relay D by TDMH. 

It will be shown later in Section IIV-AI that symmetric traffic through the relay node achieves the 

best network coding gain for MLNC. So when both buffers have packets, Algorithm \T\ schedules 

symmetric traffic; when one queue is empty certainly one-way traffic should be scheduled. By doing 

this, Algorithm Q] is able to opportunistically achieve CovH(ft TDMH , ^mlnc) by time sharing between 

TDMH and MLNC. Moreover, the attractive feature of this algorithm is that it does not require any 

arrival rate information of the two queues to achieve the system stability as stated in the following 

theorem. 

Theorem 4: Algorithm [7] stabilizes the two-way relaying system for any Poisson arrivals if the 
(bit-arrival) rate pair (Rab, Rba) £ CovH(ft TmH , 0l MLNC ) where Rab = Rab^ and Rba = Rba^- 
Proof: Without loss of generality, assume each packet has a unit length so that Rab = Rab and 
Rba = Rba- Now consider the time right after the n-th transmission. Note that each transmission 
is one of the three steps in Algorithm [Q Denote the queue length vector at time n as Q(n) = 
[lA(n) gs(n)]. It is easy to verify that Q(n) forms a non-reducible Markov chain. Note that the 
number of packets arriving at node A during a transmission time slot At is a Poisson process with 
parameter R A BAt, while the number of packets arriving at node B is also a Poisson process with 
parameter R B A<At. 

In order to show the stability of the Markov process Q(n), we define the Lyapunov function as 
follows 

V{n) 4 Emin q \{n) + Emin ql(n) + 2q A {n)q B {n), (6) 

^AB ^min ^BA ^min 

where S min = [I/Cad + ^/Cbd + l/Cmin) -1 - Now consider the case in Step 1 of Algorithm [U In 
this case we have g^(n+ 1) = ^(ra) — 1 + 5a(ti) and gs(n+ 1) = qs(n) — 1 + 8b (n), where 5^(n) 
and 5b (n) are Poisson random variables with parameters Rab /^mm and RBA/^mm, respectively. It 



10 



follows that 

E[V(n + l)\Q(n)} 



V{n) + 2- 



jAB 



R 



AB 



+2: 



J BA 



S mm V S 

Rba 



AB 



1 + 



Rba 



jab 



QA{n) 



SbA — S m i n V YiBA J ^BA 



1 + 



Rab 



q B {n)+A a , (7) 



where A a is a constant depending on Rab, Rba, S min , £ab and S#a- Next consider the case in 
Step 2 of Algorithm [TJ We have q A (n + 1) = <7a(™) — 1 + 8 A (n) and g,e(ri + 1) = 5 B {n), where 
<5a (n) and S B (n) are Poisson random variables with parameters |^ and respectively. Thus, 



E[y(n + l)|Q(n)] = V(n)+2 



Rba 



+ 



Rab 



q A (n) + A b 



(8) 



.Sab Sab ~ S m j n \Sab 

where A b is also a constant depending on Rab, Rba, S mm , Sab and S B a- 

Finally consider the case in Step 3 of Algorithm [Q We have qA(n+ 1) = <5A(n) and q B (n + 1) = 
^s(^) — 1 + 5 B {n), where 5a{ji) is a Poisson random variable with parameter |^ and 5s(n) is 
also Poisson with parameter So it follows that 



E[V{n + l)\Q{n)\ = V{n) + 2 



SbA S#a — S m j n V S^ba 



1? 



BA 



(9) 



where A c is also a constant depending on -Rab, Rba, S mm , Sab and S#a- Since (Rab, Rba) is 
within the quadrilateral formed by the origin, (£ab,0), (S 
inequalities are obvious: 



M ,a,S min ) and (0,S B a), the following 



R 



AB 



I I + ^ < and 



^ab — S min VSab / Sab ^ba — S min \Sba / S#a 



-Rba 



-Rab 



< 0. 



Therefore, we have E[V(ra + l)|Q(n)] < V(n) — 1 when qA(n) or q B (n) is sufficiently large. 
According to the Foster-Lyapunov criterion [31], Q(n) is stable. ■ 
Next, we consider the opportunistic packet scheduling algorithm for PLNC. We want to show that 
any rate pair in Theorem[3]is stabilizable. Consider a pair (Qa, Qb) £ equal to (maxQ/, maxQ b ), 
where Q f ,Q b G N+ are subject to Q f < \J1ab/£, Qb < A 2 S B a/^, and Q^AiSaba = Qf^^BAB 
1 Let (Saba, Sbab) € ^plnc be the closest point to (Saba, S B ab) with constraint QaA 2 S B ab ~ 



4 Theoretically, we should choose a positive integer pair {Qf,Qb) large enough such that (Q/^/Ai, Q&f/Aa) = (Saba,Sbab)- 
However, such a pair may not exist. So it is acceptable to select a pair with appropriate large values of Qf and Qt such that 



(Q f £/\i,Qb£/\2) » (S 



ABA : ^BAB 
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Qb^i^aba- Also, choose Q* E N + large enough such that the following inequalities are satisfied: 

A . Rab . Qa ^ n sin \ 

1 + — h — < 0, (10a) 



^ABA 


( Rba 


^BA — ^BAB 


V^BA 


^BAB 


( Rab 


Sab — SabA 


KJ^ab 



^BA Q* 

1] ^Rba + Qb < Qj 



Sab Q 

max{Q A ,Q B } - Q* < 0. (10c) 

Such Q* must exist because the region confined by (llOal) and (|10fcl) in the first quadrant is enclosed 

by the achievable rate region specified in Theorem [3] We propose the following algorithm that 

opportunistically schedules packet transmissions over the relay node by PLNC. 

Algorithm 2: Opportunistic PLNC and Scheduling 



Step 1 (qA > Qa, qs > Qb)- Qa packets from node A and Qb packets from node B are sent over 

Relay D by PLNC. 



Step 2 (q A < Qa, qB > Qb) 
Step 3 (q A > Qa, qB < Qb) 
Step 4 (q A < Qa, qB < Qb) 



min{gs, Q*} packets from node B are sent over Relay D by TDMH. 
minlg^, Q*} packets from node A are sent over Relay D by TDMH. 
qA packets from node A and qs packets from node B are sent over 
Relay D by PLNC. 



It will be also discussed in Section IIV-AI that PLNC achieves its maximum network coding gain 
when the traffic pattern over a relay node is ji = So when the two buffers respectively have 
at least Qa and Qb packets, Algorithm [2] schedules the two-way traffic such that it has a traffic 
pattern w ^j- If one °f me two buffers does not have enough packets to achieve fi w §^"> 
then one-way traffic is scheduled. However, if both of the buffers do not have enough packets then 
PLNC is adopted because it has a better throughput than TDMH for any traffic pattern as shown in 
Fig.lH 

Theorem 5: Algorithm [2] stabilizes the two-way relaying system for Poisson arrivals with the 
(bit-arrival) rate pair (Rab, Rba) within the region constructed by (|10al) - (ll0cl) where Rab = Rab^ 
and Rba = Rba^- 

Proof: The proof is omitted here since it is similar to the proof of Theorem HI ■ 
Remark 3: As shown in Theorem [5l Algorithm [2] is able to stabilize the queues at the source 
nodes with Poisson arrivals whose rate pair is within the rate region constructed by (I10al) - (|10Z?I) in the 
first quadrant, which is enclosed by IRplnc m Theorem |3] for all Q* satisfied with (llQcl) . Therefore, 
any rate pair within 3? PL nc m Theorem [3] is stabilizable since Q* can be chosen sufficiently large 
such that the region constructed by (I10al) - (ll0fcl) approaches to ft PLNC as closely as possible. However, 
this might come with some queueing delay. 
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IV. End-To-End Sum-Rate Optimization and Diversity-Multiplexing Tradeoff 

In this section, we characterize the end-to-end sum rates of the three transmission protocols with 
or without a traffic pattern constraint. Here we only consider the sum-rate optimization problem for 
a single relay since it is straightforward to extend to the multiple relay case. Next, we characterize 
the DMTs of the three transmission protocols with multiple relay nodes. 

A. End-to-End Sum-Rate Optimization 

We first consider the sum-rate optimization problems without any constraint on traffic pattern 
for TDMH, MLNC and PLNC. Then we consider how traffic pattern influences the sum rates. In 
addition, in order to compare the throughput performance we introduce the notion of network coding 
gain, and characterize under what kind of traffic pattern the network coding gain is maximized. 
Considering Theorems |T]{3] and Fig. [3l then the following result is immediate for sum rate without 
traffic pattern constraints. 

Corollary 1: (a) The maximum sum rate achieved by TDMH is max{EAB, ^ba}- 
(b) The maximum sum rate achieved by MLNC is maxjS^s, ^bb, 2£abb} if Cdb < Cda, and if 
Cdb > Cda, then it is max {T> BA , ^aa, 2£arb}- (c) The maximum sum rate achieved by PLNC is 
maxjE^s, T, BA , T, ABA + T, BA b}- 

Proof: Since the end-to-end sum rate is R AB + Rba an d all constraints in ©, © and © 
are linear, the maximum sum rate must be achieved at a corner point of the achievable rate region 
according to linear programming theory. ■ 

For the sum rate optimization problem with a traffic pattern, we have the following. 

Corollary 2: If fi is fixed, the maximum sum rates achieved by TDMH, MLNC and PLNC are: 

R* Dm = (1 + //)/ (E^ + M Ei]j) , (11) 

Kmc = ^ + v)/{C AD + vC BD . + m&x{l^}C m t n ), (12) 

R* PLNC = (l + /i)/(C7^ + /iC7^ + max{/iC7^,C7^}) (13) 
Proof: For TDMH, its achievable rate region is enclosed in the first quadrant by line + 
= 1 according to Theorem CD So replacing Rba with /j,Rab m the line equation and solving 
for R A b, CD} is obtained by (1 + p)R AB . For MLNC, we first consider the case when C DB < C DA . 
According to Theorem [2l the boundary line of the achievable rate region is |^ + = 1 when 
pi < 1. Replacing R BA with [iR AB and solving for R AB , and then we can get the desired result for 
fj, < 1. The boundary line is + |^ = 1 when yt, > 1, and hence the maximum sum rate can 
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be found by replacing Rba with (jlRab- The case Cob > Coa can be solved by the same fashion. 
For PLNC, its proof is omitted here since it is similar to the proof of MLNC. ■ 
Given the results of the sum-rate optimization, the throughput of the three transmission protocols 
can be computed. We introduce the network coding (throughput) gain p AB (in dB scale) as 

p AB ^101og 10 || (dB), (14) 

where R A and i? B are the maximum sum rates for protocol A and protocol B, respectively. Here 
if we consider TDMH as the baseline to be compared, the network coding gains of MLNC and 
PLNC are p MT = 10 log 10 (i?* LNC /i?; DMH ) and p PT = 10\og 10 (R* LNC /R* DMll ), respectively. Note that 
it is easy to check p PT > and p PM = Ppt — Pmt > 0, which show PLNC is always superior to 
TDMH and MLNC. It is also easy to verify that p MT achieves its maximum at p = 1 and p PT 
achieves its maximum at p = by Calculus. Recall that Algorithms [TJ and [2] schedule packet 
transmissions according to the packet arrival processes at the source nodes. Their main idea is to 
schedule transmissions around such optimal points. 



B. Diversity-Multiplexing Tradeoffs of Two-Way Transmission Protocols 

In this subsection we investigate the DMTs of the three transmission protocols over multiple 
relays. Diversity gain d and multiplexing gain m in [24] are redefined here in our notation as 
follows 

i a ,. loge(7) , a R{l) 
d = — hm — ■ — — and m = lim — — , 

7^00 log 7 7^00 log 7 

where e is the outage probability of information exchange in two-way relaying, R is the end-to-end 
transmission rate between source nodes, and 7 is the signal-to-noise (SNR) ratio without fading. 
Note that e and R are not defined by the traditional fashion of point-to-point transmission. They 
are defined by an end-to-end fashion because the three transmission protocols are multihop-based 
protocols, and the system we study here is aimed at the main point of information exchange so that 
it would be fairer for each transmission protocol to declare an outage in the system when either one 
source or both source nodes cannot receive the packets they desire. Thus, the outage probability of 
transmission protocol A in the two-way relaying system is defined as 



6A =P 



£a,/ [J £a,6 



(15) 



where £ A j = {XfI A j < Rab} and £ A t = {A b / Aifc < Rba} are the outage events of forward and 
backward transmission, I A j (I A; b) is the mutual information of forward (backward) transmission and 
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{Xf,Xb : Xf, Aft G [0,1], Xf + Xb = 1} are time-allocation parameters for forward and backward 
transmission, respectively. 

We first look at the DMT problem of TDMH over multiple relay nodes. Although a similar 
problem for one-way transmission has been investigated in [28] [27], our two-way results are more 
general as we can see in the following proposition. Considering Gaussian input distribution and 
the case of relay collaboration, then the forward and backward mutual information of the two-way 
relaying system in Fig. QJb) can be shown as 

-Ttdmh,/ = -^tdmh,6 = - min {Ii, I 2 } , (16) 

where I x = log (l + 7 Y.Dev AB I^a| 2 )3 and I 2 = log (l + 7 Y.Dev AB \ h Ds\ 2 ) because the for- 
ward/backward transmission first virtually passes through a SIMO channel with receive MRC and 
then through a MISO channel with transmit MRC. Note that coefficient \ means the forward/backward 
data stream needs 2 time slots. Let Dab 7^ 0- Then we have the following result. 

Proposition 1 : If all relay nodes in Dab collaborate, then TDMH achieves the following DMT: 

( 2m \ 

d = 1 - • 77TT u 7TTT1 u T ' (17) 

V mm{(l + fi)\f, (1 + l/p)X h \J 

where m G (0, min{(l + //)A/, (1 + l///)Aj,}/2). If there is no collaboration in T>ab> then TDMH 
over D* TDMH is able to achieve the DMT in (fT71) as well, where D* TDMH G Dab denotes the optimal 
relay node for bidirectional transmission. 

The proof of Proposition Q] can be found in Appendix HI The results in Proposition Q] are actually 
similar to the DMT results of the decode-and-forward protocol in [28], [27] if we set {p = 00, A& = 
1} or {fi = 0, Xf = 1}, i.e., considering the one-way transmission case. Proposition [T] is more 
general by including the time-allocation and traffic pattern influences in two-way relaying. 

The DMTs of MLNC and PLNC basically can be derived by the same way used in the proof 

of Proposition [TJ Let us first consider the MLNC and PLNC protocols with relay collaboration. 

Although the relay nodes can achieve receive MRC, they have a difficulty broadcasting to attain 

bidirectional transmit MRC at both source nodes simultaneously. So if all relay nodes broadcast at 

the same time, then the mutual information of forward and backward transmissions in this case are 
2 r (~ ~ n . 2 



Vnc,/ = g min min jX J 2 || and J M lnc,6 = ^ min |l 2 , min jX J 2 || , (18) 
hmcj = min j/i, J 2 | and /p LNC ,& = min j/i, J 2 | , (19) 



5 To facilitate the DMT analysis here, all nodes are assumed to have the same transmit power. In fact, the DMT results are nothing 
to do with the transmit powers. 
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where I x = log (l + 7] Ed g d as h DA?), h = log (l + 7I T,Dev AB h DB\ 2 ) and coefficient § is due 
to two data streams sharing three time slots. Obviously, MLNC and PLNC cannot achieve the full 
diversity gain |Dab|- Since nodes A and B loose transmit diversity from relay nodes in the broadcast 
stage, it is an suboptimal strategy to let all relay nodes broadcast the same information at the same 
time. The better idea is to find an optimal broadcast relay for MLNC and PLNC by the following 
criteria from Lemma CD 

-^mlnc = ar § max H^AiCflfl} an d -Dp LNC = arg max C DA + C DB . (20) 



In this case, h = log(l + >y\h AD * | 2 ) and I 2 = log(l + y\h BD * I ) so that I x > I u I 2 > h- Thus (1181) 
reduces to J M lnc,/ = -^mlnc,& = |min |ji,/ 2 | while (fT9l) remains unchanged. On the other hand, if 
two-way transmission is only over an optimally selected relay node, then (fT8~l) and (fl9l) become 



2 

hmcj = ^mlnc,6 = g min{log(l + 7|/i A D* LNC |),log(l + 7|^bd* lnc |)} , (21) 

2 

IplncJ = /plnc,6 = g min {log(l + 7l^D* LNC |),log(l + 7|^d* lnc |)} , (22) 
where -Dmlnc an ^ -^plnc are me optimal relay nodes determined by 

-^mlnc = ar § max ^abb and -Dp LNC = arg max Y, ABA + (23) 

D&Dab DgT) ab 

where 2'E ABB is the sum rate within 3? MLNC and T, ABA + T, BAB is the sum rate within ^ PLNC at relay 
node D E D AB . The DMTs for MLNC and PLNC with an optimally selected relay can thus be 
derived by using (|2Tj) and (l22"l) . We summarize the DMTs of MLNC and PLNC for the scenarios 
with or without relay collaboration in Proposition [2l and its proof is given in Appendix [III 

Proposition 2: If all relay nodes in D AB collaborate to receive and broadcast at the same time, 
then the following DMT is achieved by MLNC and PLNC: 

d=l ft w — ; ; x : , ■ (24) 

2min{(l + /U )A / ,(l + l//i)AJ 

where m G (0, | min{(l + /i)A/, (1 + l//x)Ab}). If all relay nodes collaborate to receive and an 

optimal relay node is selected by (|2Q|) to broadcast, MLNC and PLNC achieve the following DMT: 



3?n 

(/ = l ' D wl 1 1 " 2min{(l + / , ) A / ,(l + l/ / ,)A b }, 1 ' (25) 



Furthermore, if an optimal relay node is selected by (1231) to receive and broadcast for MLNC and 
PLNC, then the DMT in (|25l) is also achieved. 

The results in Proposition [2] are reasonable since PLNC is no longer superior to MLNC while 
broadcasting in the high SNR regime. In addition, due to relay selection diversity, using an optimally 
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selected relay to broadcast (or to receive and then broadcast) is able to achieve the full diversity 
|D|ab- The results in Propositions Q] and [2] have been presented in Fig. |4]if we set fx = 1 for the 
three transmission protocols. First look at the results of solid lines in Fig. HI They are obtained based 
on the optimal time allocation for MLNC and PLNC with two-way relaying, i.e., A/ = A& = 0.5. 
MLNC and PLNC always have a better DMT than TDMH in any transmission case. However, the 
DMTs of MLNC and PLNC may not be better than TDMH if time allocation between forward 
and backward traffic is suboptimal. For example, suppose the four time allocation parameters of 
TDMH are Ai = A 4 = 0.01 and A 2 = A 3 = 0.49 so that \ f = \ b = a01 +°- 49 = 0.5. Thus O 
becomes d = |Dab|(1 — 2m) in this case. Similarly, if the three time allocation parameters for 
MLNC and PLNC are Aj = 0.01, A 2 = 0.49 and thus A 3 = 0.5 then MLNC and PLNC have 

A / = A1+A2+2A3 = °' 34 and Xb = 1_A / = °' 66 in this case - So ® becomes d = \T> AB \(1 - 2.2m) 
and TDMH thus has a better DMT than MLNC and PLNC, as the results of dashed line shown in Fig. 
IH For two-way transmission, an ideal transmission protocol is that it can achieve full diversity gain 
I Dab I an d multiplexing gain one. The ideal DMT line in Fig. [4] can be asymptotically approached 
if there exists a transmission protocol which is able to support N source nodes to exchange their 
packets in N + 1 time slots even when N is very large. 

V. Simulation Results 

^From the analysis in Section [IVl we know the relationship between the sum rate and the 
traffic pattern parameter and the DMTs for the three transmission protocols. In this section, we 
use simulations to illustrate how the sum rate is affected by the traffic pattern in a two-way relaying 
system, and then show how outage probabilities behave as a function of SNR. We assume that all 
nodes have the same transmit power 18 dBm. The channel between any two nodes has path loss 
exponent 3.5 and is reciprocal with flat Rayleigh fading. The distance between source nodes A and 
B is 50m, and the relay nodes are randomly dropped on a 10m vertical line whose center is located 
at the middle point between nodes A and B. 

We first look at the single relay case and consider the relay node is positioned at the middle point 
of the two source nodes. The simulation result of the network coding gains for this case is presented 
in Fig. [51 From the figure, we can observe that the maximum coding gains happen at fx ~ 1 for all 
transmission protocols, and they are seriously impacted when the traffic is very asymmetric. Note 
that the maximum of p PT indeed happens at fx because is close to unity in our simulation 

setup. As expected, PLNC has not only the best network coding gain among all the protocols but 
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also the robustness against asymmetric traffic pattern. The network coding gain of the opportunistic 
MLNC protocol, p MT is also much better than what achieved by MLNC. In Fig [51 MLNC does not 
achieve a larger rate region than TDMH when > 1 or /i < 1. For example, when ji (0.2, 6.5) 
a positive network coding gain for MLNC does not exist, i.e., p MT < 0. 

Now consider how the DMTs change with different transmission protocols over multiple relay 
nodes. In all the following simulations, we set the multiplexing gain m = 1/4, fj, — 1 and A/ = 
\b = 1/2. The outage probability of each protocol is calculated by Pfmutual information < | log 7] 
with 10 8 channel realizations. First consider all relay nodes collaborate to receive and transmit (or 
broadcast) at the same time. The results shown in Fig. [6] verify the DMTs in Proposition Q] and 
(1241) in Proposition H i.e., d = \\Dab\ for TDMH and d = § for MLNC and PLNC. For example, 
the outage probability curve for TDHM with 3 relay nodes has a diversity gain about 1.5 and all 
the curves for MLNC and PLNC have the same diversity gain abut 0.6 close to |. Fig. [7] presents 
the results of selecting an optimal relay node to receive and broadcast. In the figure we can verify 
MLNC and PLNC have d = §|1>ab| which is better than d = ~\T> AB \ achieved by TDMH. The 
diversity gains of MLNC and PLNC are obviously greater than that of TDMH for any number of 
relay nodes for selection. 

VI. Conclusion 

The fundamental limits of information exchange over a two-way relaying channel with or without 
wireless network coding was investigated. We first characterized the achievable rate regions for the 
TDMH, MLNC and PLNC protocols, and found that MLNC does not always achieve a larger rate 
region than TDMH. An opportunistic time-sharing protocol between TDMH and MLNC is able 
to achieve CovH(3^ TDMH , CR MLNC ). The rate region achieved by PLNC is always larger than those 
achieved by TDMH and MLNC since PLNC can achieve the individual broadcast channel capacity. 
We then proposed two opportunistic packet scheduling algorithms for MLNC and PLNC that can 
stabilize the two-way relaying system for Poisson arrivals. The transmission performance of TDMH, 
MLNC and PLNC is also investigated in terms of the sum rates and the DMT. The maximum sum 
rates of the three protocols with or without a traffic pattern constraint were found. We showed 
that the three transmission protocols, by using the optimally selected relay node, can achieve their 
corresponding best DMTs achieved by relay collaboration. Also, we clarified that MLNC and PLNC 
may not always achieve a better DMT than TDMH due to suboptimal time allocation. 
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Appendix I 
Proof of Proposition CD 

Proof: Consider the two-way relaying system in Fig[TJa). An outage is experienced when either 
W A or Wb cannot be decoded correctly at their destination nodes. Let £ A (£ B ) denote the event that 
the relay node nodes in T>ab cannot correctly decode Wa (Wb) and E C A (£#) denote the complement 
of £ a (£#)• According to (fT5T) and using Boole's inequality, a two-way relaying system has the 
inequality of outage probability: e TDMH < e TDMH ,/ + e TDMH ,&, where e TDMH j and e TDMHj6 are the forward 
and backward outage probabilities, respectively. Since we know e TDMH ,/ = P [£tdmh,/|£a] P[£a] + 
F[E rDmJ \E c A ]¥[E c A ], e TDMH)/ = P[£ A ] + P Wh < Rab] P[£ a ] where P[£ A ] = P 



4f /i < Rab 



^ 2 

M?fe ?/za? in the following analysis, we use notation 7 * x instead of -y x in order to clearly present 
the complex expression of exponent x. Let Rab+Rba = mlog7 so that Rab = tj- log 7. Therefore, 



-TDMH,/ 



< 2-P 




min<( y ] \h DA \ 2 , ^ \ h DB\ 2 \ < 7 * df 

DdV AB 



(a) 

< 7 *(|B AJJ |d / ) 



for large 7 and m G (0, |(1 + j«)A/), where d/ = ^x^fr^y ~~ ^) anc * ( a ) foll° ws from the notation 
definition and Lemma [2] in Appendix [TTTl Similarly, e XDMH ,& < 7 * (|D,abM&) , for large 7 and m G 
(0, \{l + l/ / u)A; ) ) where d& = {^jp^jj^ — l)- Therefore, (flTT) can be obtained for large 7 and 
m G (0, |min{(l + /i)A/, (1 + l/fi)\ b }). Next, consider there is no collaboration in T> A b and an 
optimal relay is selected to receive and broadcast. In this case (fl6l) becomes 

Itdmh,/ = /iDMH.b = ^ min {log (l + i\h AD * a S) , log [l + 7l^ T * DMH s| 2 ] } > ( 26 ) 
where the optimal relay node -Dxdmh ls selected according to the following criterion: 

-D* DMH = arg min (1/C DA + 1/C DB ) = arg max 





h-DA 


2 


hDB 


2 


\h<DA 


2 


+ 


hDB 


2 



That is to choose the relay node with max{T, ab} in Dab- Moreover by substituting (|26l) into e 



TDMH,/ 



and Rab = log 7, we can obtain 



-TDMH,/ 



< 2-P 





TDMH 


2 


han* 

TDKH 


2 


/i a n* 

1 TDMH 


2 


+ 


hf>n* 

TDMH 


2 



< 7 * cL 



(b) 



< 7^(|2) AS |ci 



where (6) follows from the fact that -D^dmh e Dab is optimal and Lemmas [2] and [3] in Appendix [TTTl 
Likewise, we can get a similar result for e TDMH ,& as shown in above. So (fTTl) can be concluded. ■ 
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Appendix II 
Proof of Proposition [2] 

Proof: Using the same definitions of £^ and in the proof of Proposition [H the outage 



probability for MLNC and PLNC satisfies the inequality: e h 



< e 



NC,/ 



+ %c b, where "NC" means 



MLNC or PLNC. Here we only prove the DMT of MLNC since the DMT of PLNC can be proved by 
following the same steps. According to the proof of Proposition \T\ it is easy to show the following 
inequality for e MLNC j: 



-MLNC,/ 



< 



+ P 



A/Ii < Rab 



P 



Xjh < R 



AB 



(27) 



where P[£a] = P [§^/A < Rab] - Let R AB = Y^Mog7 and consider the first case that every relay 
node collaborates to receive and then broadcasts at the same time. Considering (fT8~l) and using 
Lemma [2l for large 7 and m G (0, 2(1 + /i)A//3), (|27l) becomes 
1 



-MLNC,/ 



< 



12) 



AB 



n ° ad 

DeD AB 



7* (\D AB \d f ) + (a fl +C7/J 7*d/ <7*d/ 



where d f = 



3m 



2(1+M)A/ 



1 , and l/cr/ a are the variances of | Y^Dev tiAD\ 2 and | 



weV A1 



h 



BD 



3m 



2(1+1/V)A 6 



1 



Then DMT in 



respectively. Similarly, we can show that e MLNC) {, < 7*d 6 where d b 
(1241) can be arrived. Now consider MLNC with optimal relay -D* LNC selected by (|20l ) to broadcast. 
Then we have the following 



Cmlnc,/ 



p 



\h, 



< 7 -kd 



(a) 



I^dsI 2 < 7 * 



(6) 

< 7* ( |Dab|<2/ 



n * 

where (a) follows from the fact that -D„ LNC I s optimal and {/i_rb} are independent, and (b) follows 
from Lemma [2] in Appendix ITTT1 Similarly, we have e MLNCfc < 7* ^|2)ab|<4) • So (1231) is obtained. 

Next, we look at the DMT of MLNC when an optimal relay node -Dmlnc ^ s selected to receive and 
broadcast. -D MLNC is determined by (T2~3T) when 7 is large. Likewise, the first step is to calculate the 
inequality of the outage probability of the forward transmission with Rab = 3^- log 7. According 



to (12Tb and (1221) . the inequality of e MLNC ,/ is obtained as follows 
"2 



P 



'-X f \og(l + j\h AD ^J 2 ) <Rab 



(c) 
< 



n p 





\h<AD 


2 


hDB 


2 


2 1 hf)A\ 


2 + 


h>DB 


2 



< 7 * 



(28) 



where (c) follows that -D MLNC is optimal and channel gains are independent, and we thus have 
Cmlnc,/ < 7* (J2)abM/) due to Lemma |3] in Appendix ITTT1 Similarly, we also can have e MLNC ,6 < 
7 * ( \ T>AB\db )■ Accordingly, (1231) can be concluded. This completes the proof. ■ 
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Appendix III 
Lemmas for DMT analysis 

We need the following definition before proceeding to prove the lemmas in this section. 

Definition: A function (7(7) : R ++ — > R ++ is said to be exponentially equal to x (i.e., (7(7) = 7*) 
if lim 7 ^oo lo ^y = x. Similar definition can be applied to other signs, such as < or >. 

Lemma 2: Let {X k , k — 1, . . . , K} be K independent exponential random variables with respec- 
tive parameter {a k , k — 1, . . . , X} and #(7) : R ++ — > If 0(7) — > as 7 — > 00 and 0(7) is 
exponentially equal to 0^, then we have the following inequality: 

K 

£**<*(7) 



P 



fe=i 



< 



7 



A0 O 



(29) 

Proof: Without loss of generality, assume that the random sequence {X k , k — 1, • • • , if} forms 
an order statistics, i.e., {Xi > X 2 > ■ ■ ■ > X K }. Thus, the event J2k=i Xk < 0(7) is equivalent to 
the intersection event of X 1 < 0(7), X 1 + X 2 < 0(7),..., J2k=i Xk < 0(7) because X k > 0, for all 
k G [!,-•• ,K\. That is, 



P 



if 



,fc=i 



P 



n(x><0(7) 
.fe=i \j=i 



< p 



A' 



,fc=i 



(a) 



A 



np 



x k < 



k 



where (a) follows from the independence between all random variables. Since all random variables 
are exponential, then we further have 



P 



X, < 



0(7) 
k 



= 1 — exp 



0( 7 ) \ W 0( 7 ) 



(30) 



where (b) follows from the fact that X k is an exponential random variable with parameter a k 
and e~ v > 1 — y, \/y G R+. By Definition and Iim 7 _ >oo log0(7)/log7 = 0^, it follows that 



P 



Ef=i ^ < « (7) < nti = t^r^ nti < 7 



rA gfcg(7) 



rA 



A6» 



The proof is complete. ■ 
Lemma 3: Let T be a given countable finite set with cardinality |T] and V be a random vector 
set whose elements are m-tuples, independent and nonnegative, i.e., V = {V;,i G N + : V, e 
R™,Vi±Vj,i ^ j}. Suppose V* G T, V t = (V tl ,V t2 , . . . ,V tm ) T G V is an exponential random 
vector with m independent entries and 7,0(7) G 0(7), {^(7)} an d {A (7)} ^ exponentially 
equal to 0^, {a^} and {/3 ioo }, respectively. Let /(V t ) and /(Vt) be respectively defined as follows: 



/(V t ) 



nr=i^+Er=iA(7)(^)' 



and /(Vt) 



nr=i ^ 



Er=iA(7)(^) ? 
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Suppose t* = argmax tGT /(V t ) and t* = argmax tea -/(V t ). If 9(j) — > as 7 — > 00, then for 
sufficient large 7 we have 

P[/(Vf.) < 0( 7 )] < 7 l T l[ e °°+ m (' 3 max-« maX )] ) 



P 



/(v t -*)<0( 7 ) 



< 7 



(31) 
(32) 



where ct max = maXi{a ioo } and /9+ ax = max i {^ oo) 0}. 

Proof: Since all random vectors in V are independent and £* = argmax te o-/(Vt), we have 

P[/(V r ) < 0( 7 )] =P[/(V t ) < 0( 7 ),Vt G T] = l[¥[f(V t ) < 0( 7 )] . (33) 



teT 



In addition, for any t E T it is easy to show that 

f(\T \ \ Y2i = i a i(l) T/ /" ^min 



i+Er=iA(7) 



m+1 



0m(7)Vt max ^ 



where V/ 



^min 



min{V t }, V w ^ max{V,}, m ( 7 ) 4 ijgjgj^ and * t 



A Vt 



m+1 



Therefore, P[/(V t ) < 6>( 7 )] < < 0~ 1 (7)^7)] < E["iP[^^ < ^(7)^(7)] • Also, 

M(7) 



< 0m(7)^(7)] 



1 — cxp 



< 



(-7; 



(34) 



'r ++ <i>m{l)^t m tm " m <Pm (7) 

where (c) follows from e -1 > 1 — x, Wx E M + and a ti is the parameter for V ti . For sufficiently large 
7, it follows that 

1 



P[/(V t )<0( 7 )]< E 



i=l 



,(7))- m ^(7)- 



(35) 



So for large 7 (1331) can be rewritten as (BTl) . For any i e T, it is easy to show that 

nr=i ^ 



/(Vt) > 



/(v t ). 



nr=i^+Er=iA(7)(^)-- 

By considering /(Vt) with constant {«j( 7 )} (so a max = 0) and the result in (I3TI) . we thus have 



/(V*) < 0(7)1 < 7^ +m ^), and P [/(V*.) < 0( 7 )1 < JI P f/( V *) < W 



P 



Therefore, (1321 ) can be arrived by (1361) . The proof is complete. 



teT 



(36) 
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Fig. 1. (a) Two-way relaying over a single relay: Time-division multihop (TDMH) protocol needs 4 time slots, but wireless network 
coding only needs 3 time slots. Rab is the end-to-end forward rate and Rba is the end-to-end backward rate, (b) Two-way relaying 
over multiple relays: All relay nodes in Dab are available to collaborate. Thus Dab becomes a big virtual relay node equipped 
with | Das | antennas, the channels from node A (Bab) to Dab (A) become a SIMO (MISO) channel. 
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Fig. 2. PHY-Layer network coding (revised from [8], [13]) 
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Fig. 3. Achievable Rate Regions for TDMH, MLNC and PLNC. Verteices (l)-(6) are (0,E S a), (T, A ba,^bab), (Sabb^abb), 
(XUb,0), (0,Ebb) and (Eaa,0), respectively. TDMH is the triangular region with vertices (1), (4), 0. MLNC is the quadrilateral 
region with vertices (1), (3), (4), in (a) and with vertices (1), (3), (6), in (b). PLNC is the quadrilateral region with vertices (1), (2), 
(4), 0. Note that MLNC, PLNC and opportunistic MLNC achieve the same rate region with vertices (1), (3), (4), if Cda = Cdb- 



d 




- m 



Fig. 4. Diversity-multiplexing tradeoffs for the three transmission protocols (|Dab| > 1 and /j, — 1). The results of solid lines are 
the case of optimal time allocation for MLNC and PLNC, i.e., Xf = Xb = 0.5. The results of dashed lines are the case of suboptimal 
time allocation of MLNC and PLNC, i.e., Xf = 0.34 and X b = 0.66. 
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1.2 




Traffic Pattern Parameter (I 



Fig. 5. Network coding gain p vs. traffic pattern parameter fi for different transmission protocols over a single relay node. 




Fig. 6. Outage probabilities of the TDMH, PLNC and MLNC protocols over multiple relay nodes with collaboration. All relay 
nodes collaborate to reach receive and transmit MRC for TDMH. For MLNC and PLNC, all relay nodes collaborate to perform 
receive MRC and then they all broadcast at the same time. 
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Fig. 7. Outage probabilities of the TDMH, PLNC and MLNC protocols without relay collaboration. An optimal relay node is 
selected to receive and transmit/broadcast for all protocols. 



